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Previous calculation of spin torques for small-amplitude magnetization dynamics around a 
uniformly magnetized state [J. Phys. Soc. Jpn. 75 (2006) 113706] is extended here to the case 
of finite-amplitude dynamics. This is achieved by introducing an "adiabatic" spin frame for 
conduction electrons, and the associated SU(2) gauge field. In particular, the Gilbert damping 
is shown to arise from the time variation of the spin-relaxation source terms in this new frame, 
giving a new physical picture of the damping. The present method will allow a "first-principle" 
derivation of spin torques without any assumptions such as rotational symmetry in spin space. 
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Magnetization dynamics driven by electric/spin cur- 
rent in submicron-scale ferromagnets has been a subject 
of active study for a decade. They include magnetization 
reversal in multilayer systems, 1-3 displacive and resonant 
motions of domain walls in thin wires, 4 ~ 14, 34 and those 
of magnetic vortices realized in disks. 16-19 These phe- 
nomena are understood as caused by spin torques, the 
effect that conduction electrons exert on magnetization 
through the microscopic s-d exchange interaction. 

There are several kinds of spin torques. For slow and 
smooth dynamics of magnetization, n, 20 as described by 
the Landau-Lifshitz-Gilbert equation, 

n = joH eS xn + aonxn + t' eh (1) 

the following torques are known; 
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They include the celebrated spin-transfer torque (first 
term; v s being the "spin-transfer velocity" ), 21-25 the 
so-called '/3-term' (second term), 12-14 ' 26-29 the Gilbert 
damping (third term), and the renormalization of spin 
(fourth term). 30 These torques are "adiabatic" in the 
sense that they are first order in space/ time deriva- 
tive. 31,32 In particular, the a- and /3-terms, which are 
classified as "dissipative torques" 29 and arise from spin 
relaxation of conduction electrons, 12 has been a current 
controversy. 11,14,25 ~ 29 ' 35 This is due to the phenomeno- 
logical nature of the argument, and thus their micro- 
scopic understanding and derivation is highly desired. 

There are several works on the microscopic determi- 
nation of a and / g. 26 - 29 ^ 36 -38 p or ex ample, in ref. 27 (to 
be referred to as I), we have carried out a fully micro- 
scopic calculation based on the linear response of con- 
duction electrons to small-amplitude magnetization fluc- 
tuations around a uniformly magnetized state. However, 
application of the results to finite-amplitude dynamics 
requires an assumption of rotational symmetry on the 
spin torque, as exploited in writing the functional form 
of eq.(2). 
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In this paper, we present a formulation which is not 
restricted to the small-amplitude dynamics, and thus al- 
lows a "first-principle" derivation of the spin torque with- 
out any assumptions such as rotational symmetry. This 
is achieved by working with an "adiabatic" spin frame, 
where an inhomogeneous (and time-dependent) magne- 
tization is transformed to a homogeneous (and static) 
one with a price of the appearance of SU(2) gauge field. 
This method based on gauge fields has been used in, e.g., 
spin-fluctuation theory, 39 quantum tunneling and resis- 
tivity of domain walls, 40 and spin-transfer torque. 21,24 
Extension of this method to dissipative torques (a- and 
/3-terms) is the aim of the present work. 

We study the same model as I, namely, the s-d model 
consisting of conducting s electrons and localized d spins. 
The total Lagrangian is given by L to t = Ls + L c \ — H s d, 
where Ls is the Lagrangian for d spins, and 

d h 2 - 



dr c T 
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H sd = —M 



J drn{x)-{c)erc) a 
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(4) 



Here c' — c^(x) — (cj(x), c|(x)) is the electron creation 
operator at x = (r,t), and cr is the vector of Pauli spin 
matrices. We adopt a continuum description for d spins, 
with a fixed magnitude S and varying direction n. 

The s electrons are treated as a free electron gas in 
three dimensions subject to impurity potential Vj mp — 
u J^i S{i — Ri) + u s J2j &j — R'j); the second term 
represents quenched magnetic impurities, which is aimed 
at introducing spin relaxation process. The averaging 
over impurity spin direction is taken as Sf = and 







Iq2 X.xat/3 



(5) 



Namely, we consider isotropic spin scatterings, or sim- 
ply, "isotropic impurities" , to focus on the essence 
of the gauge-field formulation. More general case of 
"anisotropic impurities" [i.e., (Sf) 2 = (S*f) 2 ^ JSff 
as considered in I ] will be commented on at the end. 

To treat finite-amplitude dynamics of magnetization, 
we work with a local/instantaneous spin frame (called 
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"adiabatic frame" below) for s electrons whose spin 
quantization axis is taken to be the local/instantaneous 
d-spin direction, n. The electron spinor a(x) in the new 
frame is related to the original spinor c(x) as c(x) = 
U(x)a(x), where U is a 2 x 2 unitary matrix satisfy- 
ing a(n-cr)c = a^a z a. We take U(x) = m{x)-cr with 
m = (sin(0/2) cos </>, sin(#/2) sin0, cos(0/2)) based on 
the parametrization n = (sin 9 cos (f>, sin 6 sin <fr, cos 0). 

Since d^c = W(d^ + Ud^a = W(d M + iAJa, there 
arises from L c \ a coupling Hi = hj dr [Ao ■a- + A i -j si ] + 
0(A i ) of a-electrons with an SU(2) gauge field 

= -%rf(dJJ)=A°o a = A^<T, (6) 



An 



Here A^ = m x d^m is a measure 



in the new frame 
of temporal (/x = 0) or spatial (// = 1,2,3) variation of 
magnetization, and couples to spin, cr(a^) = (a^era)^, and 

spin current, j s i — {H/2mi)[a'a , dia], in the new frame. 
The spin-torque density from H s d is given by 

t eX {x) = Mn x (er(x)) nc = M1Z (z x {a{x)) ne ), (7) 

where <r(x) = (c'<tc) x , and (■ • -) nc means averaging in 
nonequilibrium states with time-dependent n or steady 
current. We have introduced a 3 x 3 orthogonal matrix 

K al3 = 2m a m 13 - 6 af3 , (8) 

representing the same rotation as U but in a three- 
dimensional vector space. Note that IZz — n, c^crc = 
ft(aW), and K 2 = 1. 

Since the gauge field A^ contains a space/time deriva- 
tive of magnetization, one may naturally formulate a gra- 
dient expansion in terms of A^ to calculate, e.g., the 
torque (or spin polarization). In particular, the adiabatic 
torques are obtained as the first-order terms in A^: 

2 r 

+ b^(z x A M ) 



M 



(9) 



Here a± = a — z (z-cr) and A^ = A^ — z (z-A^) are the 
respective transverse components, 41 and the sums over 
fi = 0, 1, 2, 3 are understood. From the identities, 

KA^-^nx^n), Tl(z x A x ) = -c^ra, (10) 

together with eq.(7), we see that eq.(9) leads to the 
torque density t c \ = a^h + (a ■ V) n + bo (n x n) + 
n x (6 -V) n, which is related to t' cl of eq.(2) via t c \ — 
(HS/a 3 ) tgj, with a 3 being the volume per d spin. 

The coefficients a M and o M are expressed by response 
functions which are wavevector (q) and frequency (lu) 
dependent, and the relation between (<x±) nc and A^ is 
nonlocal in general. For adiabatic torques, however, it is 
sufficient to estimate them at q = and u> = 0, since the 
gauge fields already contain space/time derivative. The 
torque then becomes a local function of n like eq.(2). 

The following calculations are done with Green's func- 
tions, Gfc CT (z) = (z — h 2 k 2 /2m + £fct + i7o-sgn(Imz)) _1 
where 7 CT = h/ (2r CT ) = 7m;u 2 iv + Tm a u 2 S 2 mp {2v s + v a )/3 
is the damping rate, and e-p a = £f + &M . Here n x (n s ) is 
the concentration of normal (magnetic) impurities, and 
v a = m k~Fcr /2tt 2 H 2 (with k-p a — y^Smifv/S) is the den- 
sity of states at £f<t- We calculate the torques in the 
lowest nontrivial order in 7< t /£fo' and j^/M, as in I. 42 



Gilbert Damping: We first study the torques asso- 
ciated with the temporal variation of n. Assuming n to 
be uniform, n — n(t), we focus on the \x = terms in 
eq.(9). The coefficients ao and bo are related to the trans- 
verse spin susceptibility, x°l '> as Jii~ a ^°L + bo£ af3 ) = 



foXjfil = 0,uj = 0), whereas we obtain (Fig.l) 



27 
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(11) 



where p s = n^ — n^ (with n a = fc Fcr /67r ) is the s-electron 
spin density, S'f = 8 af} - 5 az 5P z , e xy = -e yx = 1, and 

£ xx = £ yy = o Thig leadg tQ ss = psCl 3/ 2 an d a = 0. 

The Gilbert damping does not arise. One may question 
the result (11), but it seems robust since it assures the 
fact that, when the direction of n changes slightly, the 
equilibrium s-electron spin density follows it with the 
same magnitude p s . Thus the gauge-field method has an 
apparent difficulty with Gilbert damping. 

This difficulty can be resolved if we note that the im- 
purity spins, which are static (quenched) in the original 
frame, become time-dependent in the adiabatic frame. 
Namely, the spin part of Vimp is expressed as 



Sj-c'crc = Sj(t)-a' era 
where Sj(t) = ^(t)" 1 ^ = K(t)S 3 , or 



(12) 



(13) 



is the impurity spin in the adiabatic frame, which is 
time-dependent. While the s-electron spin tends to fol- 
low the instantaneous d-spin direction n(t), it is at the 
same time pinned by the quenched impurity spins Sj. 
These two competing effects are represented by the time 
dependence of Sj(t) in the adiabatic frame. The above 
observation (and difficulty) is based on the neglect of this 
time dependence. 

To study the effect of time variation of Sf (t) , we cal- 
culate (ct" (t)) perturbatively with respect to Sj(t). The 
nontrivial contribution comes up at the second order 
(since Sf — 0). Using the extension of the linear response 
theory to the second-order response, 44 we obtain 



(52(*))x 



dt' 



dt" X ^Jt-t',i/-t") 



x SP(t')Si(t"), (14) 

where jg^t - t',t' - t") = -(2ft 2 )- 1 9(t - t')9(t' - 

t")([[aZ(t),aP(t')},^(t")}) + (f - t",0 «- j) vs the 
retarded three-spin correlation function, here playing the 



a 1 




a" 



+ 




Fig. 1. Transverse spin susceptibility, Xjf> U P to 0(-/a-), con- 
tributing to the spin torque related to the time variation of mag- 
netization. The thick (thin) solid line represents an electron line 
carrying Matsubara frequency ie n + iu\ (ien). The dotted line 
with a cross represents scattering by static (quenched) impuri- 
ties, either non-magnetic or magnetic. 
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role of nonlinear response function. 45 From eqs.(5) and 
(13), we have 

I , .„ , .. „ , >- (15) 



The correlation function x^p 1 (t — t',t' — t") will decay as 
a function of time separations, with a time scale of the 
spin relaxation time r s . For lot s <C 1, where w is a typical 
frequency of magnetization dynamics, we may make a 
short-time approximation for the time-dependence of n 
(and TV), and retain the lowest nontrivial order in the 
expansion, TZ(t') = K(t) + (*' - t)TZ(t) + ■■■. This leads 
to 

TZ(t')n(t") =l-(t'- t")K(t)li(t) + • • • , (16) 

where we have used !Z(t)lZ(t) = —!Z(t)lZ(t). The identity 

[Jl(t)li(t)} af) = 2e a ^Al (17) 

shows that we have found a contribution of the gauge 
field, which was not expected at the initial stage of the 
formulation when we introduced the gauge field, eq.(6). 
We rewrite eq.(14) in terms of Fourier components, 



(^i(w))nc = n s u s 



, , °° du>' R 



2tt 



and expand (tjj,tu') with respect to u> and ui' as 



(18) 



(19) 



up to the first order, where A a p 7 = Xa/3 7 (0j 0); an d -B a/ 3 7 
and C a j3j are the coefficients. Back to real time, we have 



d-~r. 



(52(i)) M - n s u 2 s \ A a01 S?{t)S~<{t) + B afiy —SP(t)S^(t) 



= -n s u 2 s S? mp {A aW + 2C aM s^ x A*{t) + ■■■} 

(20) 

Hence it suffices to calculate the w'-linear coefficient, 
CVr> of X^ 7 (w,w'). 

The retarded function Xa(j 7 (w,w') can be obtained 
from the Matsubara function 

Xaf}7 (iuA,iu'x) = l f d(T-r") f d{r'- T ") 
1 Jo Jo 



x e 



iu> x (T-T')+iu>' x (T' -t") 



(T T [a Q (r)a /3 (r')a 7 (r")]) 



(21) 



by the analytic continuation, iu>\ — > u> + 2irj and ilo' x — > 
to' + irj, where r\ is a positive infinitesimal. It is suffi- 
cient to consider the lowest-order process (with respect 
to the spin scattering), whose diagrammatic expression 
is shown for (cr" (iLJ\)) ae in Fig. 2. The corresponding 
Xa0j(iu\,iuj' x ) is given by 



+ iuJ\)cr 



n k.k' 



x G k '(ie n + i^)cr 7 Gfe(ie„)J + (iuj' x — > iw A - iu^), 

(22) 

where G k — diag(Gfc|, Gki)- To leading order in 7 CT , we 
have 

= ^ E tr [° aG l° P K' -G$)<fGi] -((3 <-> 7 ), 

k.k' 

(23) 

and thus C aP . ( = -{■K/2M)vl(8'f5~i z - S^S?*), where 
v + = i/f + . Therefore, we obtain 

2tt 
~3M' 

leading to the Gilbert damping with 



7 n s u 2 s Sf m y + (zxA^(t)), (24) 



7T 9 ^2 2 ^ 

" = 3^^^^+ x — , 



(25) 



in agreement with the result of I for isotropic impurities. 

The above calculation provides us a new picture of 
Gilbert damping. In the adiabatic frame, impurity spins 
are dynamically fluctuating and disturb the s-electron 
spins. Because of this disturbance, s electrons admit a 
low-energy spectral weight in the transverse spin channel 
(at wavevector q = 0), which causes Gilbert damping. 

Current-Induced Torques: We next calculate the 
spin torques induced by electric/spin currents in the pres- 
ence of spatial variation of n = n(r). These torques cor- 
respond to terms in (<x_i_) n e proportional to the spatial 
components of the gauge field, Af~ (and z x Af). The 
texture n (and hence Ai) is assumed to be static here. 

We produce a steady current state by applying a d.c. 
electric field E, and first calculate the linear response of 






Fig. 2. Contribution to the Gilbert damping in the lowest order 
in 7 CT . The wavy line represents scattering from impurity spins, 
which are dynamically fluctuating in the adiabatic frame. The 
running frequencies are indicated. 




Fig. 3. The Ei- and A -linear coefficient, K a , of the transverse 
spin polarization, (ct" (q)) nc , in the presence of current flow (Ei) 
and magnetization texture (A?). The vertices with Vi and a^Vj 

are associated with Ei and , respectively. Other graphical 
meanings are the same as Fig.l. 
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cr" to E. Using Kubo formula, 44 one can derive 

,,„, ., ,. Kf (q, u> + iO) — Kf(q, 0) ^ .„„. 
{at (q)) nc = lim -L^l -> Ei , (26) 



K?{q,iu x ) 



dre^ (T T al(q,T)Ji). (27) 



Here J, = -eJ2k v * c l c k = - e Efc^4°fe ~ 
(eh/m) J2k, q ' ^q'.i a I+q' cr/3a fc' witn w i = hki/m, is the 
total electric (or charge) current, and the average (• • •) 
is taken in the thermal equilibrium state determined by 
L c \ — H s d- We extract the gauge field in the first order, 
and write as 



Kf(q,iw x ) 



(28) 



Diagrammatic expressions for Kf? are shown in Fig. 3. 
The calculation is rather straightforward, yielding 



a 



2e Js 



b 



2tt 
3M 



■n s u. 



2 g2 



'2e 3s 



(29) 



Here j s — a s E = J-j- — j± is the spin current, with <j s = 
(e 2 /m)(n|T|— njrj,) being the "spin conductivity" . 47 The 
result (29) also agrees with I for isotropic impurities. 

Before closing, let us comment on the case of 
anisotropic impurities. For the Gilbert damping, we can 
reproduce the result of I if we take the anisotropy axis 
to be n. 46 As for the /3-term, we could only reproduce 
the first term of the previous result [eq.(25) or (26) of I], 
and not the second term. There seems to be a key con- 
tribution to the torque from 14 mp , but we do not succeed 
in reproducing the result of I at present. While this puz- 
zle needs to be clarified, we would like to stress that the 
essence of the gauge-field method is already contained in 
the isotropic case as presented in this letter. 

In summary, we have developed a theoretical scheme 
for the adiabatic spin torques in terms of gauge fields 
associated with the adiabatic spin frame of conduction 
electrons. We have found an unexpected contribution of 
the gauge field from the source terms of spin relaxation, 
which is essential to give the Gilbert damping. The mes- 
sage here is that the gauge-field method does work also 
for dissipative torques (a- and /3-terms), if we are care- 
ful enough and faithful to the microscopic model. This 
method will be readily applicable to other type of spin- 
relaxation mechanisms such as the ones originating from 
various type of spin-orbit coupling. 
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